Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 17 January 2011 



(MN MgX style file v2.2) 



Magnet ot her mal and magnetorotational instabilities in hot 
accretion flows 

De-Fu Bu 1 *, Feng Yuan 1 ! and James M. Stone 2 

C ' 1 • 1 Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, 
£S1 ' Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, China 
^ 2 Department of Astrophysical Sciences, Peyton Hall, Ivy Lane, Princeton, NJ 08544> USA 

hj ; 

. 17 January 2011 

PQ ■ ABSTRACT 

In a hot, dilute, magnetized accretion flow, the electron mean-free path can be much 
greater than the Larmor radius, thus thermal conduction is anisotropic and along 

■ magnetic field lines. In this case, if the temperature decreases outward, the flow may 

be subject to a buoyancy instability (the magnetothermal instability, or MTI). The 

Q . MTI amplifies the magnetic field, and aligns field lines with the radial direction. If the 

S— j ' accretion flow is differentially rotating, the magnetorotational instability (MM) may 

r/i \ also be present. Using two-dimensional, time-dependent magnetohydrodynamic simu- 

. lations, we investigate the interaction between these two instabilities. We use global 

simulations that span over two orders of magnitude in radius, centered on the region 
. around the Bondi radius where the infall time of gas is longer than the growth time 

J> ■ of both the MTI and MM. Significant amplification of the magnetic field is produced 

t-H ' by both instabilities, although we find that the MTI primarily amplifies the radial 

C<~) . component, and the MM primarily the toroidal component, of the field, respectively. 

CO ■ Most importantly, we find that if the MTI can amplify the magnetic energy by a factor 

| F t , and the MM by a factor F r , then when the MTI and MM are both present, the 

— . magnetic energy can be amplified by a factor of F t ■ F r . We therefore conclude that 

t— I ' amplification of the magnetic energy by the MTI and MM operates independently We 

f^*) \ also find that the MTI contributes to the transport of angular momentum, because 

. radial motions induced by the MTI increase the Maxwell (by amplifying the magnetic 

' field) and Reynolds stresses. Finally, we find that thermal conduction decreases the 

• i-H , slope of the radial temperature profile. The increased temperature near the Bondi 

1 radius decreases the mass accretion rate. 

■ Key words: accretion, accretion discs - conduction - magnctohydrodynamics: MHD 

- ISM: jets and outflow - black hole physics 



1 INTRODUCTION 

Hot, diffuse accretion flows are likely present in low- 
luminosity active galactic nuclei (AGN) and the quiescent 
and hard states of black hole X-ray binaries (Narayan 2005; 
Yuan 2007; Narayan & McClintock 2008). When the mass 
accretion rate is very low, the electron mean-free path will 
be very large, so the accreting plasma is nearly collisionless. 
For example, from Chandra observations of the accretion 
flow at the Galactic Center, Sgr A*, the electron mean-free 
path is estimated to be 0.02 — 1.3 times the Bondi radius 
(~ 10 — 10 6 Schwarzschild radii). In such cases, thermal 
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conduction can have significant influence on the dynamics 
(Quataert 2004; Johnson & Quataert 2007), resulting in the 
transport of thermal energy from the inner (hotter) to the 
outer (cooler) regions. If the energy flux carried by ther- 
mal conduction is substantial, the temperature of the gas 
in the outer regions can be increased above the virial tem- 
perature. Thus, gas in the outer regions is able to escape 
from the gravitational potential of the central black hole and 
form outflows, significantly decreasing the mass accretion 
rate (Tanaka & Menou 2006; Johnson & Quataert 2007). 

Time-dependent hydrodynamical simulations of hot ac- 
cretion flows show that they are convectively unstable, no 
matter whether the flow is radiatively inefficient (Igumen- 
shchev & Abramowicz 1999; Stone, Pringle & Begelman 
1999; Igumenshchev & Abramowicz 2000) or efficient (Yuan 
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& Bu 2010). This is consistent with the prediction of one- 
dimensional analytical models of advection-dominated ac- 
cretion flows (ADAFs; Narayan & Yi 1994; 1995). Convec- 
tive instability occurs because the entropy of the flow in- 
creases in the direction of gravity. However, in a magnetized 
weakly collisional accretion flow, in which the electron mean- 
free path is much greater than the Larmor radius, thermal 
conduction is anisotropic and primarily along magnetic field 
lines. In this case, Balbus (2000; 2001) found that the sta- 
bility criterion for convection is fundamentally altered: con- 
vection occurs when the temperature - not the entropy - 
increases in the direction of gravity. Convective instabilities 
in this regime are referred to as the magnetothermal insta- 
bility (MTI). 

In a weakly magnetized, non-rotating atmosphere, lo- 
cal simulations have demonstrated that the MTI can am- 
plify the magnetic field and lead to a substantial convective 
heat flux (Parrish & Stone 2005; 2007). Sharma, Quataert 
& Stone (2008; hereafter SQS08) investigated the effects of 
the MTI on accretion flows in two-dimensions using global 
simulations in spherical polar coordinates. Their simulations 
focused on the region around the Bondi radius, where the 
MTI growth time is shorter than gas infall time. They found 
the MTI can amplify the magnetic field and align field lines 
with the direction of the temperature gradient (radially), 
consistent with the results obtained by previous local simu- 
lations. The effects of the MTI on the intracluster medium 
(ICM) have also been investigated (Parrish, Stone & Lemas- 
ter 2008). In that study also, the MTI drives the magnetic 
field lines to be radial, leading to a conductive flux which is 
an appreciable fraction of the classical Spitzer value. Such 
a high rates of thermal conduction make the slope of the 
temperature profile flatter. 

The studies of the MTI described above all focus on a 
non-rotating accretion flow. It is well-known that in a weakly 
magnetized, differentially rotating hot accretion flow, the 
magnetorotational instability (MRI; Balbus & Hawley 1991; 
1998) will be present, and MHD turbulence driven by the 
MRI is thought to be the most promising mechanism of an- 
gular momentum transport. In principle the MTI and MRI 
should coexist in such flows. Therefore, an interesting ques- 
tion for investigation is whether the MTI and MRI interact 
with one another, i.e., does one instability enhance or sup- 
press the other, or do the two instabilities operate entirely 
independently. For example, it is of interest to ask whether 
convective motions driven by the MTI can transport angular 
momentum, and similarly whether turbulent motions driven 
by the MRI can transport heat. 

Motivated by the question above, we present the results 
of two-dimensional magnetohydrodynamical (MHD) simu- 
lation of hot accretion flows with anisotropic thermal con- 
duction and rotation. We investigate flows which are slowly 
rotating (not rotationally supported) at the Bondi radius, 
where the growth time of both the MTI and MRI is shorter 
than the infall time. This paper is organized as follows: we 
introduce our numerical methods and initial conditions in 
Section 2, most of our results are presented in Section 3, 
while Section 4 is devoted to a summary and discussion. 



2 NUMERICAL SIMULATION 

We adopt two-dimensional, spherical coordinates (r, 9) and 
assume axisymmetry d/d<j> = 0. We use the Zeus-2D code 
(Stone & Norman 1992a, 1992b) to solve the MHD equations 



with anisotropic thermal conduction: 

| + pV.v = 0, (1) 

P^ = -Vp- P V^+i-(VxB) xB, (2) 

^ = Vx(vxB), (3) 

P^ = - P V.v-V.Q, (4) 

Q = - x b(b ■ V)T. (5) 



Here p is the mass density, v is the velocity, p is the pressure, 
ip — —GM/r is the gravitational potential (G is the gravita- 
tional constant and M is the mass of the central black hole), 
B is the magnetic field, e = p/(7 — 1) is the internal energy 
(7 is the adiabatic index), Q is the heat flux along magnetic 
field lines, % ls the thermal diffusivity, b is the unit vector 
along magnetic field and T is the temperature of the gas. 

For convenience, we adopt \ = np/T and assume 
that k = a c (GMr) 1 ^ 2 which has the dimension of a 
diffusion coefficient (cm 2 s _1 ). Anisotropic thermal con- 
duction is implemented using the method based on limiters 
(Sharma & Hammett 2007) to guarantee that heat always 
flows from hot to cold regions. 

We initialize our simulations as follows. The radial ve- 
locity, internal energy, and density of the accretion flow are 
adopted from the hydrodynamic Bondi solution. For the 
other two components of the velocity (vg and v$), we set 
vg = and v<f, = lo/r with lo(= const.) (a parameter of our 
simulations) the specific angular momentum of the accretion 
flow. The initial magnetic field is uniform and perpendicular 
to the equatorial plane, (i.e. B = B z , where z is the vertical 
cylindrical coordinate). We have run four models. The pa- 
rameters used in these models are listed in Table 1. In models 
LI and L2, we set lo to be 0.1 times the Keplerian angular 
momentum at the inner boundary. Thus, in models LI and 
L2, the centrifugal barrier is inside our inner boundary, and 
rotationally supported region is never formed in these two 
models. In models HI and H2, however, we set lo to be 10 
times the Keplerian angular momentum at the inner bound- 
ary. Thus, in models HI and H2, the centrifugal barrier is 
outside our inner boundary, so that a rotationally supported 
disk can form, which is subject to the MRI. We set j3 (the 
ratio of gas pressure to magnetic pressure) to be f3 — 10 12 
at the outer boundary in our initial conditions. We assume 
that the adiabatic index 7 = 1.5. We use the Bondi radius, 
tb = GM/c%> (coo is the sound speed at infinity and we 
set c 2 ^ = 10 _6 c 2 where c is the speed of light) to normalize 
all of the length scales in our paper. Times reported in this 
paper are all in the unit of the inverse rotation frequency at 
the outer radius of the domain, fi~Jt = [r^ ut /GAl] ll/2 . 

All of our calculations span a domain bounded by an 
inner boundary at Ti n = 0.05r_e and an outer boundary 
at 20rs. As in SQS08, we adopt a numerical resolution of 
60 x 44. The radial grid is logarithmically spaced. A uniform 
polar grid extends from 6 = to 7r. At the poles, we use 
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Table 1. Parameters for all of our models 



Models 


ll 


al 


description 


LI 


O.lhn 





no MTI, no MRI 


L2 


0.1l in 


0.2 


only MTI 


111 


Wkn 





only MRI 


112 


10hn 


0.2 


with MRI, with MTI 



t Iq is the initial angular momentum of the flow. li n is the 
Keplerian angular momentum at the inner boundary. 
X = K P/T is the thermal diffusivity and ft = a c {GMr) 1 ^ 2 . 



axisymmetric boundary conditions. In the radial direction, 
we use inflow/outflow boundary conditions at both the inner 
and outer boundaries. For the temperature at the inner and 
outer boundaries, we use outflow boundary conditions. 



3 RESULTS 

3.1 Magnetic field amplification 

A weakly magnetized, conductive, differentially rotating ac- 
cretion flow is subject to magnetothermal and magnetoro- 
tational instabilities when the condition, 



1 dp d\nT dQ, 2 Q 
p dr dr dlnr 



(6) 



is satisfied, where Q. is the angular velocity of the accretion 
flow. Provided that the magnetic field is weak, the condition 
for the existence of MTI alone is 



1 dp dlnT 



< 0. 



(7) 



p dr dr 

This condition is always satisfied in the initial conditions 
of our calculations, thus as long as the MTI growth time 
is smaller than the gas infall time, the MTI will grow. The 
condition for the existence of the MRI is 

^<°- < 8 > 
dmr 

This condition is again always satisfied in the initial condi- 
tions of our calculations. Defining 

1 dp d\nT 



Imti = 



JMRI 



p dr dr 

dn 2 



dlnr 



(9) 
(10) 



the MTI and MRI growth timescales are then t mti = 7^ T7 
and tuiRi = Jmri' respectively. The wavelength of the max- 
imum growth rate for MTI and MRI are 



\MTI = 2irV a /"fMTI 

and 



(11) 



XMRI = iKVahMRI, (12) 

respectively. Here v a (= \f B 2 /farp) is the Alfven speed. 

We investigate the interplay between the MTI and MRI 
by examining the magnetic field amplification by the two in- 
stabilities. Note that in addition to the MTI and MRI, the 
magnetic field can also be amplified by the geometric con- 
vergence of the accretion flow in spherical geometry and flux 
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Figure 1. Evolution of the volume averaged magnetic energy in 
both the r and (f> components of the magnetic field normalized 
to the initial magnetic energy in a shell from r" out /40 to r out in 
model Ll. In this model, magnetic field is amplified only by flux 
freezing. Since vg ~ 0, there is almost no amplification of the 6 
component of the field, so it is not shown. The growth of each 

{{B 2 r,e,<j>) v ~ B r,B,4,,o)> where 



component is defined as SB 2 , 



,0 



Qv denotes volume average and Bq is the initial field strength. 



freezing. First, we investigate the magnetic field amplifica- 
tion by flux freezing by carrying out model Ll. We note 
that in models Ll and L2, the centrifugal barrier is inside 
our inner boundary, thus the timescale for the MRI to grow 
is much longer than the gas infall timescale. Furthermore, 
in model Ll we set the coefficient of the thermal diffusivity 
a c — 0, i.e., there is no MTI in model Ll. Therefore, the 
magnetic field amplification in model Ll can only be due to 
flux freezing. The results below show that the effect of flux 
freezing on magnetic field amplification is very small, allow- 
ing us to ignore this effect in the analysis of models L2, HI 
and H2. 

In model L2, we set the coefficient of the thermal dif- 
fusivity a c = 0.2 so that it is unstable to the MTI. In this 
model, the magnetic field is amplified only by the MTI. In 
models HI and H2, the centrifugal barrier is outside our in- 
ner boundary. The timescale for the growth of the MRI is 
much shorter than the gas infall time, so that the MRI will 
affect the flow in both of these two models. In model HI, 
we set the coefficient of the thermal diffusivity a c = 0, so 
that the MTI is suppressed and the magnetic field is ampli- 
fied only by the MRI. In model H2, we set the coefficient 
of the thermal diffusivity a c = 0.2, so that both the MTI 
and MRI coexist, and the magnetic field is amplified by the 
nonlinear evolution of both. The properties of the four mod- 
els are summarized in Table 1. The models demonstrate the 
interplay between the MTI and MRI. We have also tried 
various other values of a c , and find that changing a c by a 
factor of a few does not affect our results, consistent with 
the conclusion in SQS08. 

Figures 1-4 show the change of magnetic energy as a 
function of time in these four models. The energy is averaged 
from r = r out /40 to r out . The solid, dotted, and dashed lines 
correspond to the energy in the r, 9, and (f> components of 
the magnetic field, respectively. 
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Figure 1 shows the magnetic energy evolution in model 
LI. The field is smoothly amplified by flux freezing. At the 
end of the simulation, the field is only amplified by a factor 
of 2. Because vg ~ 0, there is almost no amplification of the 
9 component of the field. From this model, we can see that 
the effect of flux freezing on magnetic field amplification is 
very small. 

We now investigate the magnetic field amplification by 
the MT1 (model L2). The only difference between models LI 
and L2 is that in the latter we include thermal conduction 
so that the MTI is present. Fig. 2 shows the evolution of the 
magnetic energy in model L2, it shows significantly more 
amplification of the magnetic field in comparison to model 
LI. This amplification is mainly due to motions driven by 
the MTI. Initially, the field is rapidly amplified by the MTI. 
After time t = 6, the magnetic field ceases to grow expo- 
nentially, and the growth rate becomes small, for the follow- 
ing reason. The growth rate of the fastest growing mode of 
the MTI is given by 7 ~ \[l/p][dp/dr][dlnT/dr]\ 1/2 b e where 
bg = Bg/~B. The growth rate decreases with time because 
of the decrease of bg with time. In this model, we use pa- 
rameters F r ,t, Fg tt and F^ : t to represent the amplification 
factor of the r, 9 and </> components of the magnetic field, 
respectively. We find 

F r ,t = Blj/Bl ~ 100, (13) 

Fg,t = Bl f /Bl ~ 10, (14) 

F*,t = Blj/Bl ~ 100, (15) 

where B r j, Bgj and B$j correspond to the r, 9 and cj) com- 
ponents of the magnetic field at the end of the simulation, 
respectively. The radial component of the magnetic field is 
amplified by the MTI significantly. The MTI saturates by 
aligning field lines with the temperature gradient (radial 
direction), which is consistent with the results of previous 
work (Parrish & Stone 2005; 2007; SQS08). SQS08 found 
that the magnetic field amplification by MTI is quasi-linear. 
If the magnetic field is weak initially, MTI can amplify the 
field by a factor of ~ 10 — 30, independent of the initial field 
strength. Our calculation confirms their result. 

We investigate field amplification by the MRI in model 
HI. In this model, we turn off thermal conduction so the 
MTI does not exist, and the magnetic field amplification 
must be due to motions driven by the MRI. Figure 3 
shows the time evolution of the magnetic energy. Initially 
(t < 16.5), the field is very weak, and at our resolution the 
fastest growing mode of the MRI is not resolved, so that field 
amplification occurs primarily by flux freezing. The initial 
(t < 1) rapid amplification of the azimuthal component of 
the magnetic field is due to shear. At t = 16.5, the magnetic 
field strength is amplified by flux freezing to a high enough 
value that the MRI is resolved, and the field undergoes sub- 
stantial and rapid amplification by the MRI. To illustrate 
this point, we define a parameter Re = Xmri/Ax, where 
Aa; is the grid size, which measure how well the MRI is re- 
solved on our numerical grid. Figure 5 shows contours of Re 
at t = 16.5. From this figure, it is clear that unstable wave- 
lengths of the MRI are well resolved at that time. The MRI 
amplifies the field rapidly, and shortly after t = 16.5 the 
MRI saturates and the field begins to grow algebraically. 
Saturation of the MRI can be understood as follows. The 
condition for MRI to operate is k 2 v 2 + dQ, 2 /d\nr < 0. After 
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Figure 2. Change of volume averaged magnetic energy in the r, 
8 and <j> components of the magnetic field normalized to the initial 
magnetic energy in a shell from r ou t/40 to r out as a function of 
time in model L2. In this model, the magnetic field is amplified by 
the MTI. The growth of each component of the magnetic field is 
defined as SB 2 0ij> = (<B 2 r ,e,«p) v - B 2 g 4> Q ), where ()y denotes 
volume average and Bg is the initial field strength. 
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Figure 3. Change of volume averaged magnetic energy in the r, 
8 and <j> components of the magnetic field normalized to the initial 
magnetic energy in a shell from r ou t/40 to r Q ut as a function of 
time in Model HI. In this model, magnetic field is amplified by 
the MRI. The growth of each component of the magnetic field is 
defined as &B 2 6(j> = ((B 2 r> g }<p ) v - B 2 g 4> Q ), where ()y denotes 
volume average and Bois the initial field strength. 

t = 18, we find k 2 v 2 + dQ, 2 /dlnr ~ in the region where 
the MRI growth wave length can be resolved, so the growth 
rate of MRI decreases significantly. Using parameters F r , r , 
Fg -r , and F$ fT to represent the magnetic field amplification 
factor of the r, 9 and </> components of the magnetic field, 
respectively, we find 

F r , r = Blj/Bl ~ 10, (16) 
F e>r = Blj/Bl - 10, (17) 
F<t,, r = Blj/Bl ~ 10 4 , (18) 
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Figure 4. Change of volume averaged magnetic energy in the r, 
9 and <j> components of the magnetic field normalized to the initial 
magnetic energy in a shell from r out /40 to r out as a function of 
time in Model H2. In this model, magnetic field is amplified by 
both the MTI and MRI. The growth of each component of the 



magnetic field is defined as SB 2 , g ± = {{B 2 r g l f > ) v 
where ()y denotes volume average and Bq is the initial field 
strength. 




becomes algebraic. Using parameters F r , Fg and to rep- 
resent the amplification factor of the r, 9 and <j> components 
of the magnetic energy, respectively, we find 



F r = B r j/B 

Fe = B 2 j/Bq 



10 , 
100, 



F = Blj/Bl ~ 10 6 



(19) 
(20) 
(21) 



where B r j, Bgj and B^j correspond to the r, 6 and cj> com- 
ponents of the magnetic field at the end of the simulation, 
respectively. 

Comparing models L2, HI and H2, we find, 



F r ~ F r .t ■ F r .r, 

Fe ~ Fg t • Fg r , 



F A 



(22) 
(23) 
(24) 



Let us first consider the amplification of the radial compo- 
nent of the magnetic field. If the MTI alone amplifies the 
radial component of the magnetic energy by a factor of F T> t 
(model L2), and the MRI alone by a factor of F r . r (model 
HI), then when the MTI and MRI coexist the radial com- 
ponent of the magnetic field can by amplified by a factor 
of F r ,t • F r , r (model H2). The amplification factors for the 
other two components of the magnetic energy behave simi- 
larly. We therefore conclude that in the nonlinear regime the 
MTI and MRI amplify the field independently, most likely 
because the MTI is buoyancy instability driven by radial 
gradients in the temperature, while the MRI is a dynam- 
ical instability driven by the radial gradient of the orbital 
frequency. Provided that the magnetic field is weak and the 
plasma is dilute, thermal conduction from hotter to cooler 
regions will make the plasma buoyantly unstable, regardless 
of whether the flow also contains turbulent eddies driven 
by the MRI. On the other hand, provided that the mag- 
netic field is weak and differential rotation exists, the MRI 
can operate regardless of the presence of convective eddies 
driven by the MTI. 
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Figure 5. Contour of Re = Xmri I Ax at time t = 16.5, where 
Ax is the grid size. It is clear that the wavelength of the fastest 
growing mode of the MRI is well resolved at time t = 16.5. 



where B r j, Bgj and B^j correspond to the r, 9 and cj> com- 
ponents of the magnetic field at the end of the simulation, 
respectively. The toroidal component is most significantly 
amplified as expected. 

In model H2, we turn on thermal conduction so both 
the MTI and MRI coexist. Figure 4 shows the time evolu- 
tion of the magnetic energy. Initially (t < 3), there is only 
the MTI. This is because initially the field is too weak for 
unstable wavelengths of the MRI to be resolved. At t = 3 
the magnetic field has been amplified by the MTI enough 
that unstable modes of the MRI are resolved, and the MRI 
begins to drive motions. After t = 3, the field begins to be 
amplified exponentially by both the MTI and MRI, until at 
t ~ 20 both instabilities saturate and the growth of the field 



3.2 The role of MTI in angular momentum 
transfer 

The MRI is now believed to be the most promising mecha- 
nism of angular momentum transport in fully ionized accre- 
tion flows. It is interesting to investigate whether the MTI 
can contribute to angular momentum transfer in hot, diffuse 
accretion flows. 

We first compare the angular momentum profiles in 
models LI and L2 at the end of both calculations. The flow 
in model LI is laminar, while there is (only) MTI in model 
L2. Angular momentum profiles at the end of the simula- 
tions are shown in Figure 6. It is clear that the angular 
momentum distribution of model LI is almost the same as 
the initial distribution. This is because the flow in model 
LI is laminar, and the magnetic field is very weak: therefore 
both the Reynolds and magnetic stresses are very small and 
there is almost no angular momentum transfer. It is also 
cleat that in model L2, angular momentum is transferred 
from the inner to the outer regions. In order to investigate 
the mechanism for the angular momentum transport in this 
case, we plot the Maxwell and Reynolds stresses in Figure 7. 
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Figure 6. Radial profiles of specific angular momentum for mod- 
els LI (dotted line) and L2 (solid line). The dashed line shows the 
initial angular momentum distribution. 



Figure 7. Radial profiles of stress. The dotted and solid lines 
correspond to the Maxwell stress in models LI and L2, respec- 
tively. The triangle and square are Reynolds stress in models LI 
and L2, respectively. 



The dotted and solid lines correspond to the Maxwell stress 
(—BrB^/in) in models LI and L2, respectively. Because of 
magnetic field amplification by the MTI, the Maxwell stress 
in model L2 is much larger than that in model LI. The tri- 
angles and squares denote the Reynolds stress in models LI 
and L2 respectively. Due to the presence of motions driven 
by the MTI in model L2, the Reynolds stress is three orders 
of magnitude higher than in model LI. Because the mag- 
netic field is still very weak in model L2, the Maxwell stress 
is much smaller than the Reynolds stress. We conclude the 
angular momentum transport in model L2 is dominated by 
the Reynolds stress. We note that although the MTI induced 
Reynolds stress does transfer angular momentum, the effect 
is quite low, so that there is only a small change of the an- 
gular momentum profile compared to the initial conditions. 

Next we compare the angular momentum distribution 
in models HI and H2 at late times in Figure 8. Both profiles 
are quite similar. Recall that both models are unstable to the 
MRI, while only H2 is unstable to the MTI. The similarity 
in the angular momentum profiles indicates that the role of 
the MTI in angular momentum transport is less important 
compared to the MRI. 

3.3 Effects of MTI on the structure of the 
accretion flow 

Figure 9 shows the time-averaged density (upper-left panel), 
temperature (upper-right) and mass accretion rate (lower 
panel) for models HI and H2. In the upper panels, the solid 
and dashed lines correspond to models HI and H2, respec- 
tively. In the lower panel, the solid, dotted and dashed lines 
correspond to the net accretion rate, inflow rate and out- 
flow rate, respectively. The black and red lines correspond 
to models HI and H2, respectively. As shown in the upper- 
right panel, the radial profile of temperature in model H2 
(with MTI) is flatter than that in model HI. This is likely 
because thermal conduction transports energy from the in- 
ner to the outer regions of the flow. From the lower panel, 
we can see that the inflow rate decreases due to the effects 
of thermal conduction. This is consistent with Johnson & 
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Figure 8. Radial profiles of specific angular momentum for mod- 
els HI (solid line) and H2 (dashed line). 

Quataert (2007). The decrease is caused by the increased 
heating of the outer regions by thermal conduction. 



4 SUMMARY AND DISCUSSION 

In extremely low accretion rate systems, the plasma is very 
dilute, the collisional mean-free path of electrons is much 
greater than their Larmor radius; thus thermal conduction is 
anisotropic and along magnetic field lines. Anisotropic ther- 
mal conduction can fundamentally alter the Schwarzschild 
convective instability criterion. Convection sets in when 
temperature (as opposed to entropy) increases in the di- 
rection of gravity. Convective instability in this regime is re- 
ferred to as the magnetothermal instability (MTI). The MTI 
can amplify the magnetic field and align field lines with the 
temperature gradient (i.e., the radial direction). If the flow 
has angular momentum, the MRI may also exist. We have 
investigated the possible interaction between the MTI and 
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Figure 9. Radial profiles of time-averaged variables in models 
HI and H2. Upper-left and upper-right panels correspond to the 
radial profiles of density and temperature, respectively. In the 
upper panels, solid and dashed lines correspond to models HI 
and H2, respectively. In the lower panel, the solid, dotted and 
dashed lines correspond to the net accretion rate, inflow rate and 
outflow rate, respectively. The black and red lines correspond to 
models HI and H2, respectively. 

MRI by performing two-dimensional magnetohydrodynam- 
ical (MHD) simulations of nonradiative rotating accretion 
flows with anisotropic thermal conduction. 

We have presented the results from four models, with 
focus on the amplification of magnetic energy by the com- 
bination of MTI and MRI. As a control case, we computed 
a model (LI) in which neither the MTI nor the MRI was 
present, so that the magnetic field is amplified only by flux 
freezing and geometrical compression in the flow. We find 
that field amplification via flux freezing in this case is very 
small. Next, we investigated the effects of the MTI on mag- 
netic field amplification in model L2. We find that the MTI 
can amplify magnetic field significantly, which is consistent 
with the result obtained in previous work. Finally we consid- 
ered the evolution of rotating flows that are unstable to the 
MRI. In model HI (no thermal conduction), only the MRI 
is present. In model H2, both the MTI and MRI coexist. 
By comparing the structure and evolution of these differ- 
ent models, we find that if the MTI alone can amplify the 
magnetic energy by a factor of Ft, and the MRI alone by a 
factor of F r , then when the MTI and MRI are both present, 
the magnetic energy can be amplified by a factor of Ft ■ F r . 
Thus, we conclude that the MTI and MRI operate indepen- 
dently, and that there is little feedback on the amplification 
factor in the nonlinear regime. Because the MTI and MRI 
are driven by different properties of the flow (temperature 
gradient in the case of the former, and angular velocity gra- 
dient in the case of the latter), the presence of one may not 
strongly affect the presence of the others. Thus, anisotropic 
thermal conduction can make the plasma buoyantly unsta- 
ble, regardless of whether or not turbulent eddies driven by 
the MRI are present. Similarly, an unstable rotation profile 
can result in the MRI, regardless of whether or not there are 
convective eddies induced by the MTI. 

In addition to the interplay between the MTI and MRI, 
we have also investigated the effects of MTI on angular mo- 



mentum transport. We find that the MTI enhances angular 
momentum transport simply because it can amplify mag- 
netic field (thus enhancing the Maxwell stress) and induce 
radial motions (thus enhancing the Reynolds stress). But 
we note that the effects of the MTI on angular momentum 
transfer are small. Finally, we find that thermal conduction 
does make the slope of the temperature smaller due to the 
outward transport of energy, and the increase of the tem- 
perature in the outer regions decreases the mass accretion 
rate. 

Lastly, we would like to mention that in a weakly col- 
lisional accretion flow, the ion mean-free path can be much 
greater than its Larmor radius; thus the pressure tensor is 
anisotropic. In this case, the growth rate of the MRI can in- 
crease dramatically at small wavenumbers compared to the 
MRI in ideal MHD (Quataert, Dorland & Hammett 2002; 
Sharma, Hammett & Quataert 2003). In this regime, the 
viscous stress tensor is anisotropic, and Balbus (2004; see 
also Islam & Balbus 2005) has shown that when anisotropic 
viscosity is included, the flow is subject to the magnetovis- 
cous instability (MVI). The MVI may also amplify magnetic 
field significantly. An interesting project for the future is to 
investigate the effects of the MVI on hot accretion flows. 
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